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Set Books (Paperback Editions) 


D. L. Kreider, R. G. Kuller, D. R. Ostberg and F. W. Perkins, An Intro- 
duction to Linear Analysis (Addison-Wesley, 1966). 


E. D. Nering, Linear Algebra and Matrix Theory (John Wiley, 1970). 


It is essential to have these books; the course is based on them and will 
not make sense without them. 


Conventions 


Before working through this correspondence text make sure you have read 
A Guide to the Linear Mathematics Course. Of the typographical conven- 
tions given in the Guide the following are the most important. 


The set books are referred to as: 


K for An Introduction to Linear Analysis 
N for Linear Algebra and Matrix Theory 


All starred items in the summaries are examinable. 


References to the Open University Mathematics Foundation Course Units 
(The Open University Press, 1971) take the form Unit M100 3, Operations 
and Morphisms. 


4.0 INTRODUCTION 


One of the principal applications of linear mathematics is to linear dif- 
ferential equations. You have seen already in the Foundation Course 
(Units M100 24 and 31, Differential Equations I and IT) how linear differential 
equations can arise in the theories of population growth and mechanical 
vibrations. Another application is to electric circuits; it is the subject of the 
television component of this unit. There are many other applications, 
mostly in physical science and the technology based on it; you will meet 
some of these later in the course. 


In general, differential equations are harder to solve than algebraic equa- 
tions because the solutions are functions rather than numbers. However, 
just as an algebraic equation is simpler to solve when it is linear, so too 
a differential equation is easier to solve when it is linear; because it is then 
a linear problem (in the sense defined on page N63), and the solution set 
then has a particularly simple algebraic structure. For example, in the 
differential equation 


Sf" +f=0 (domain of f is R$) 


which we used in the Foundation Course (Unit M100 31, Differential 
Equations IT) to model a mass-spring system, the solution set consists of 
all functions obtainable by giving real values to A and B in the formula 


Jf: tı Acost+Bsint (te R3). 


ZW 
+10 
! die 
3! 0 
8: 
s -5 
i -10 


In other words, the solution set is a vector space, consisting of all linear 
combinations of the functions cos and sin; in the notation of Unit 1, 
Vector Spaces, it is the set spanned by cos and sin, i.e. 


(cos, sin). 


Why is the solution set of the differential equation a vector space? The 
reason is that the equation can be written in the form 


L(f)-0 


where L is a linear transformation of some vector space of differentiable 
functions; the solution set of such an equation is the kernel of L, and we 
know already (page N31) that the kernel of any linear transformation is a 
vector space. In the present case the linear transformation in question is 
defined by 


N=" f 


with a suitable domain. 
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In this unit we will investigate the ideas sketched above and generalize 
them to yield more general linear differential equations (i.e. differential 
equations that are linear problems). In particular, we shall see how to 
reduce the problem of solving any linear differential equation of first 
order to the evaluation of integrals. 


In order to be able to complete the solutions to these and other problems 
involving integration, you will need to use some techniques which are 
based on the rules of integration by substitution and integration by parts 
(see Unit M100 13, Integration II). 


It is not the purpose of this course to teach these techniques, but the main 
results are collected in the supplementary handbook, Techniques of Inte- 
gration (referred to as TI) to which you should refer when necessary. 


41 LINEAR DIFFERENTIAL OPERATORS 
4.1.0 Introduction 


1n solving the differential equation 
J"+f=0 

we are ]ed to consider the linear transformation L defined by 
L(f)=f“+f  (fesome domain) 


Such a linear transformation, mapping any function in its domain to a 
linear combination of itself and one or more of its derived functions is 
called a linear differential operator. We shall give a more precise definition 
later, but before we can do so, it is necessary to decide what set of functions, 
to use as the domain for such an operator. 


4.1.4 The D Operator 


The simplest linear differential operator is the differentiation operator 
Difv——f' — (fe some domain of D). 


We know that D is a linear transformation from previous work (Unit 2, 
sub-section 2.1.5, Exercise 3) but what about its domain? The obvious 
choice for the domain of D is the set of all real functions (functions whose 
domain and codomain are R or subsets of R) that are differentiable every- 
where in their domain. (This means that f'(x) exists for all x in the domain 
of f.) We shall call such functions real differentiable functions. The obvious 
choice is not always the best, however, and here two departures from the 
obvious are desirable. First, in order to solve equations of the form 


Df — (given function in image set of D) 


we would like to make use of the Fundamental Theorem of Calculus. This 
theorem (as stated in Unit M100 13, Integration II) tells us that differentia- 
tion can be reversed by integration provided the function to be integrated 
is continuous. So, to make sure that the fundamental theorem will apply, 
we would like all the functions in the image set to be continuous. We can 
do this if we re-define the domain of D to consist of the real differentiable 
functions f, such that f” is continuous. The graph of such a function, f, is 
not only continuous but is also smooth: the slope of the tangent varies 
continuously. Such functions are said to be continuously differentiable. 
(Perhaps you are wondering why we require f to be continuously differenti- 
able: surely the derived function of a differentiable function is necessarily 
continuous? Unfortunately this is not so; there are “ pathological" dif- 
ferentiable functions whose derived functions are not continuous. An 
example is given in the Appendix.) 


integrate 
Continuously 
differentiable Continuous 
functions functions 
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A second departure from the obvious is to specify the type of subset we 
will admit as the domain of our function f: so far we have implied, rather 
imprecisely, that any subset of R would do. The Fundamental Theorem of 
Calculus works for real continuously differentiable functions, provided the 
domain has no gaps. A subset of R without gaps is called an interval; for 
example, R and R* are intervals. You have already met notations such as 


[a, b] for the interval (x: a < x < b) 


in the Foundation Course. We use the notation C(I) for the set of all real 
functions continuous in some interval J; the particular case C[a, b] was 
used in Unit 1, Vector Spaces as an example of a vector space (see page 
K3). We also use the notation C'(J) for the set of all functions continuously 
differentiable in the interval J; this set is a suitable domain for the linear 
differential operator D. 


ct) c(t) 


Example 


We have introduced a few words and notations in this section: to help 
them sink in we look at an example. 


Let f be the modulus function 
> 
fix—| x zd (xe R) 
Clearly f is continuous: so f e C(R). But f is not differentiable because it 
has no derivative at 0. 
Suppose we restrict the domain of f to [1, 3], say. Then fe C[I, 3] and 
fi:x-— 1 (xe[1,3]) 


So f“ is continuous and, hence, fis continuously differentiable in [1, 3]; i.e. 
fe 1, 3]. 


As further examples, the sine and cosine functions belong to C(R) and 
CHR). 
Exercises 


1. Let be the function defined by 
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f(x) 


f(x) t 


1 (x <0) 
x-—(1+x (0<x<1l) 
0 (x» 1). 
Which of the following statements are true? 
(a) fe C(R) (b /e C[—1,1] 
(g) fe C'[- 1, 1] (d) fe C'D, 4]. 


2. Give an example of a function not already given that is continuous but 
not continuously differentiable in the interval [— 1, 1]. 
3. Is it true or false that C'[a, b] c G[a, b]? 


4. Verify that C! [a, b] is a real vector space. (Hint Look at the subspace 
criterion at the bottom of page K12. We discussed this in Unit 1.) 


Solutions 
l. (a) FALSE (b) TRUE (c) FALSE 
(d TRUE 
2. One example is the function 


p EID) een 


f(x) f(x) 


1 1 


gap in domain 


The function f is continuous, since its graph has no breaks, 
but it is not differentiable at 0, since 


fQ-5-f0 [r1 if A>0 
h -[ 0 if a 


so that iB t) O does not exist (see Unit M100 7, 
hoo 


p. 30). 
(If you were unable to find an example yourself, try and find 
another one, now. It is a good habit to adopt for this type of 
exercise.) 

3. True; for every differentiable function is continuous. 
Proof If the function f is differentiable at a then the limit 
in the formula 


f'@)= is 0-70) 


exists. 
To prove that f is continuous at a, we have to prove that 


lim f(x) 

x78 
exists and is equal to f(a). 
We can write 

lim f(x) = lim f(a +h) 
whence ™ 


lim f(x) — f(a) = lim f(a + h) — f(a) 
h^0 


xa 


=limh x [erento 

h^-0 
qo 
=0 x f'(a) 
=0 


4. The solution is given in Example 3 on page K13. Be sure to 
read it. 


* This step follows from the fact that if lim f(x) and lim g(x) exist then lim (f(x) g(x) exists 


and is equal to their product. 


4.1.2 Powers of the D Operator 


For differential equations of the second and higher orders we shall need, 
in addition to D, operators giving the second and higher derived functions 
of a function. The simplest of these is the operator mapping a function f to 
its second derived function f"; we saw in the Foundation Course that this 
operator is the composition of D with itself, i.e. 


Do D(f) = D(D(f)) =f". 


It is usually abbreviated to D?. 


To define D? properly we should specify its domain, as we did for D. 
Once again, if we are to invert the operator D? using the fundamental 
theorem twice, we shall want the functions we integrate to be continuous 
—that is, we want both f" and f” to be continuous. 


continuous continuous 


In fact, it is sufficient to require f" to be continuous; for then ° is con- 
tinuously differentiable, and any differentiable function is continuous (see 
Exercise 3 of sub-section 4.1.1). Any function f whose second derived 
function f" is continuous is said to be twice continuously differentiable. We 
denote the set of twice continuously differentiable functions in some 


2, 
c*() cel) ca) 


interval I by C^(I). Thus the domain of D? may be taken as ©), with 7 
any specified interval of R. (We proved that D? is a linear transformation 
in Exercise 3 of sub-section 5.3 of Unit 2.) 


Example 
The function 
fix-—x? (xeR 


belongs to C?(I) for any interval I. 


f'(x) 


f(x) f(x) 
D 
—— 
r T—— —r > —p — > 
x x 
fix — x2 
fix r+ 2x 


On the other hand, consider the function 


. x? x20 
sx—[ 5 22 (xe R) 
g(x) 
xex? 
x 
xi -x? 
graph of g 


This function is differentiable everywhere: the only point requiring in- 
vestigation is 0. The graph would seem to indicate that there is a derivative 


11 


at 0, but one could also verify this formally by considering 


lim £C + ñ) — £O 
Bm nha 


A0 h 


for h positive and ^ negative: in either case the limit is zero. Also g is 
continuous: therefore, g e G1( R). 


g(x) 


Xi -2x x> 2x 


graph of g’ 


Clearly g" is defined except at 0, So, for instance, we have 
g¢C(R), g¢ C[-1, 1], ge CU, 3]. 


In a similar way we can define linear transformations D?, D*, etc.; their 
domains can be taken as C?(I), C*(I), etc., for any I where C*(I) is the set 
of all functions three times continuously differentiable in the interval J, 
and so on.* 


For example, f e C"(R) means that the nth derived function of fis continu- 
ous at all points of the real line. 


It is an immediate consequence of the definition that each of these spaces, 
-.. (I), C(I), CHI), CW), is a subset (in fact a subspace) of its pre- 
decessor; in symbols we write 


CHa CD > CN > CU) >- 


Exercises 
l. (i) Is the function x ——— x? + x3, (x e R) an element of 
(a) C([0,1, (b C7f0,1], (o) @25[0,1]? 
(i) Is the function x ——— |x|x?, (x e R) an element of 
(a) CXR), (b) CR) (9 eo» 


2. Given an example of a function belonging to C'(R) for every n. 

3. Givean example of a function belonging to C' [ — 1, 1] butnot G?[-1, 1]. 
(Hint Can you use Solution 2 of sub-section 4.1.1?) 

4. Give an example of a function belonging to C'[—1, 1] but not 
to C^*![—1, 1]. 


Solutions 
l. (i) (a) Yes, (b Yes, (c) Yes. 
Gü) An alternative definition of the function is 


—[ 5 a (xe R) 


* Strictly, the domain of a function is an integral part of the function, so that D" can only have 


onc domain. But we assume that D" is given a specific domain C"(1) in any context, and we 
do not invent a new symbol at cach occurrence. 
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and this is now similar to the example we discussed 
above. 


The derived function is 


3? x20 


* [x x«0 


) (xe R) 
and we know, from the example, that this latter function 


belongs to C'(R). Therefore the original function belongs 
to C?(R). 


The answers are, therefore, 
(a) Yes, (b) Yes, (c) No. 


2. The zero function is the simplest example, though perhaps 
too trivial for some tastes. x ———+ x? + x° isalso an example. 
Among the many others are constant functions, powers 
(x —— x” with ne Z*), polynomial functions, the expo- 
nential function, the sine and cosine functions. The space of 
functions with this property is denoted by C*(I). We may 
loosely say that such functions are “infinitely ” differentiable 
on J. i 

3. If fis the required function, the condition to be satisfied is 
equivalent to the statement that f' is continuous but not 
continuously differentiable. By solution 2 of sub-section 4.1.1, 
one such function is given by 


fe-[ iazo Cete 


if x«0 


Since this f" is continuous we can obtain f by integration; we 
obtain 


m fy? i 
7-[r-E p Exo) @eDL1D 


where k is any constant. The figure depicts the case k —0. 


f(x) f) f'(x) 


There are many other solutions to this exercise; e.g. 
2 

_ {x 1>x>0 

To= jo ucc 


4. We can generalise the method used in Solution 3: f™, the 
nth derived function of f, is to be continuous but not con- 
tinuously differentiable; so we may take f(x) to be x if x 
> 0 and 0 if x < 0. Integrating z times we obtain 
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f= (n+ D! (e[-1, 1) 


where & is any constant. 
Once again there are many alternative answers. 


This example shows that each of the inclusion relations in the 
chain mentioned earlier 


eu) > C) > CW 3 :::, 


is proper; i.e. for each n € Z* there are functions in C'(I) 
which are not in C"*!(I). 


4.1.3 Polynomials in D 


The linear transformation we considered at the beginning of this unit, 
Lif——f"-^f  (fedomain of L) 
can be expressed as the sum of two simpler linear transformations 


y— 
and 

fF 
with the same domain. (The sum of linear transformations was discussed 
in Unit 2, sub-section 1.5. See also Section 2-2 of K.) The first of these is 
the linear transformation we have been calling D?, and the second is the 


identity transformation; thus it is consistent with our previous notation to 
write 


L=D +1 


This gives a very convenient notation for linear differential operators; for 
example, the transformation 


LSS mof 


(with domain C?(I) for some interval I) is denoted by 3D? — D + 1. 
Differential operators of this type are important because differential 
equations of the form 


Af= g, 


with A an operator of this type, arise quite frequently in applications and 
can always be solved in principle. In the following reading passage this 
class of operators, the so-called constant-coefficient linear differential 
operators, is defined formally. 


READ Section 2-3 on pages K48-K53 omitting Example 5. 


Much of this reading passage is revision of work we studied previously in 
N. 


Notes 


o line 5, et seq. page K51. We shall not be interested in the concept of nilpotence 
ere, 


Gi) line ~8, page K51. “. . . factorization of polynomials depends only...” 
This bit is not strictly true. If you examine in detail all the properties of ordinary 
arithmetic that you use in order to write 

(x + 2)(x + 3)x — 1) = x° + 4x! + x — 6, 


you will find that there are others than those mentioned by K. Fortunately, the 
are all properties which linear transformations possess. cd 
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(iii) line —6, page K51. Examples of functions in Cla, b] are 
x—! (x e [a, b] 
X — —5 exp x (x e [a, bD) 


(iv) line 5, page K52. K does not distinguish between functions and images. 
In the notation of the Foundation Course, which does make this distinction, we 
would either regard y as a function throughout, and write 


(D? + D—2Dy =y" + y — 2y 


or else regard it as a symbol for the image of x under a function denoted by some 
other letter, say Y, in which case the equation would have to be written 


(D° + D —2)Y)G) = Y) + Y'6) — 2YG) 
or 


d'y | dy 
2 " gp ee A 
(D? + D — 2) Y)G) dat dx 2y 


When we come to manipulate specific functions, the notation used in K has its 


advantages. The authors' explanation of their notation is given at the end of this 
reading passage. 


Exercise 


Exercise 7, page K54; parts (a), (b), (c) and (d). 


Solution 
The answers are given on page K728. The method is a straight- 
forward application of the definition; e.g. (using the K notation 
for functions) 


(D? — 1)2e* = D'Qe*) — 1 - 2e* = 2e* — 2e* = 0 


For (c) and (d) there are two alternative methods either using the 
bracketed operators in succession, or multiplying their product out 
first 


(D + 1)(D — 2)e* = (D + 1)(De* — 2e*) = (D + 1)(-e*) 
= -e — e” = — 2e* 
or 
(D + 1)(D — 2)e* = (D? — D — 2)e* 
= (D'e* — De* — 2e*) 
= e* — e — 2e* = —2e* 
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4.1.4 Linear Differential Operators in General 


To be able to deal with as wide a class of linear differential equations as 
possible, it is useful to consider linear differential operators more general 
than the constant-coefficient type defined in the previous sub-section. This 
generalization is given in the next reading passage. 


READ Section 3-1 from page K86 as far as “. . . the operator itself.” on line 
6 of page K88. 


Notes 


(i) line 4, page K86. By x D? + D + x is meant the operator which maps any 
real function f e C?(7) to the real function 


x — —95 xf (x) + f“(x) + xf (x) (xe1) 


(See also Example 5, page K52.) 

Gi) Definition (3-1), page K86. As an example, xD? + D + x is a linear dif- 
ferential operator, of order 2. The choice of C"(/) for the domain ensures that D, 
D?, ..., D" all produce continuous functions, and multiplying these continuous 
functions by the continuous functions ao, 4;,..., a,, gives further continuous 
functions, so that all the images under L, obtained by adding these continuous 
functions, are themselves continuous. However the operator y ——— xy” + yy’ 
+ y!? cannot be expressed in the form of Equation (3-1) and is not linear. 
(ii) line —4, page K86. “‘a,(x) is not identically zero on J" means “‘a,(x) is not 
equal to zero for all x € /" or in other words *' a, is not the zero function on I”. 
(iv) lines I and 3, page K88. The notation used in K for intervals is defined in 
the footnotes on pages K3 and K88. Here is a full specification. 


For a, b € R, we have 


[a, b] = (x:a < x < b):a and b both lie in the interval 

(a, b) = (x:a < x < b): neither a nor b lies in the interval 

[a, b) = (x:a < x < b):a, but not b, lies in the interval 

(a, b] ={x:a <x < b): b, but not a, lies in the interval. 
These intervals are usually referred to as being “closed”, “open” and “ half- 
open" respectively. Unfortunately the symbol (a, b) has two meanings: the 
interval defined above, and the ordered pair whose two elements are a and b. 


(That is why we used the notation Ja, b[ for open intervals in the Foundation 
Course.) 


Where relevant, we write + eo (“ plus infinity”) for b and — co for a. This is just 
a notational trick; there are not points actually called + co or — eo in R. Thus, 
when the symbols co and — co are used, they go with a round bracket rather than 
a square bracket, to indicate that they are not themselves included in the intervals 
they define. We have the following possibilities: 


[a, +0) = (x: a < x) 
(a, + co) = (x: a < x) 
(— co, b] = (x: x < b) 
(— œ, b) ={x:x <b} 
(—%, + co) = R 
(v) line 2, page K88. The coefficients of Dz, D, etc., can be any continuous 


functions. All that is required for linearity is that the operator is equivalent to 
the form of Equation (3-1). The notation In denotes the natural logarithm function. 


Exercises 


l. Exercise 1(c) and I(d) on page K88. If you require further practice try 
Exercises 2, 3 on K88, 89. 


2. State the order of the linear differential operators in Exercise 2 page 
K88, on the interval (0, 1]. 


16 
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Solutions 


1. The answers are given on page K731. 
In (d) you must apply the bracketed operators one at a time 


to the function. The fact that 
(D+1(D-x» #0? +(1-x)D-x 


is taken up in the next sub-section. 
2. (a) 2ifa#0;lifa=0,b#0;0ifa=b=0. 
(b) 2. (c) 2. 


4.1.5 Products of Linear Differential Operators 


We saw earlier (top of page K52) that constant-coefficient linear differen- 
tial operators can be multiplied just like ordinary polynomials, so that, for 
example 

(D+2XD-1)=(D-XD+2)= D? + D —2. 


For general linear operators, however, multiplication does not follow this 
simple law, and in fact it is not even commutative. An example illustrating 
the procedure for multiplying such operators is given in the next reading 
passage. 


READ the remainder of Section 3-1 on page K88. 


It is crucial to remember that an operator L is not a disembodied expression 
but has a meaning according to its effect on a typical element of its domain. 
Whenever in doubt about re-expressing the form of L, write out L(y) in 
full. 


Exercises 


1. Verify the non-commutativity of the multiplication of two general 
linear differential operators by carrying out calculations similar to the 
one in the reading passage, for: 


(i) the product (xD + 2)(3D + 1) and 
(ii) the product of the same two operators in the reverse order. 


2. Exercise 4, parts (a), (b), (c) on page K89. 


Solutions 
1. (i) A calculation similar to the one on page K88 gives 
(xD + 28D Dy = (xD + 2)3y' + y) 
= 3xy" + xy’ + 6y' + 2y 
= (3xD? + (x + 6)D + 2)y 
(ii) On the other hand 
(3D + 1)(xD + 2y = (3D + Go! + 2y) 
= xy" + 3y + 6y 
+ xy’ + 2y 
= (3xD? + (x + 9)D + 2)y 
thus (xD + 2)8D + 1) =3xD? + (x + 6)D + 2 
but (3D + D(xD + 2) 2 3xD' + (x +9)D +2 
2. (a) (DP+1(D-1) =D?-D?+D-1 
Being constant-coefficient lincar differential operators, 


(D? + 1) and (D — 1) can be multiplied like ordinary 
polynomials. 


(b) Following the method on page K88, we apply the dif- 
ferential operators to an arbitrary function y in an appro- 
priate vector space (€7(J) will do here, though it is not 
in general necessary to specify this). 

(xD(D — x)y = xD(y' — xy) 
= xy" — x?y' — xy 
ie. xD(D—-x)-xD'—xiD-x 

(g) (xD? + D/P f= (xD? + Df" +f) 

-xDGf?  f* f?) + SOS +f 

= xf) + fO 2/0) + Cf + 2f") 

= x(xf + 4/07) + 2f" 

Therefore (xD? + D)?f = x? D* + 4x D? + 2D’. 

In writing out operators with non-constant coefficients, it is essential to 
write the coefficients (x, — x?, — x, in part (b) above) to the Jeft of the D 
symbols (just as in Equation (3-1) on page K86). This is to keep them out 
of the line of fire of the D operator, which acts on the functions written to 
the right of it. 


4.1.6 Summary of Section 4.1 


In this section we defined the terms 


continuously differentiable (page K13) soe 
interval (page K3) — 
linear differential operator (page K86) see 
order (page K86) e.. 


We introduced the notation 


[05] (page K86) 
Cla, b] (page K3) 

D?, D" (page K50) 
@(D (page K86) 
(a, b) (page C16) 
(a, co) (page C16) 
(—co,b) (page C16) 
(00,00) (page C16) 


Techniques 

1. Use (aD + b)(cD + d) = acD? + (ad + bc)D + bd when a, b, c and t. 
d are constant functions. 

2. When a, b, c and d are not constant functions, apply the operator in | ` `. 


by applying each factor in turn. 
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42 LINEAR DIFFERENTIAL EQUATIONS 
4.2.0 Introduction 


You have already met some linear differential equations in the Foundation 
Course; for example the equation 


f'*f-0  (fiR,* — R) 


mentioned earlier in this unit, which models the vibration of a mass-and- 
spring system. Using the notation we have developed in this unit the equa- 
tion can be written in terms of the linear differential operator D? + 1, in 
the form : 


(D? -Df-0 
where fe C?(R,*). 


Thus the problem of solving the linear differential equation is equivalent 
to finding the kernel of the linear transformation D? + 1 (the inverse 
image of the element 0 in the codomain). Such an equation is called a 
linear problem by N (page N63); K’s terminology is operator equation (page 
K80). 


This point of view makes it possible to exploit the general theory of linear 
transformations in solving such differential equations. You have already 
met all the main ideas of this section in the Foundation Course (Units 
M100 24 and 31, Differential Equations I and II); all we do here is to 
revise these ideas, define the technical terms more precisely, and extend 
them to the general linear differential equation. 


4.2.1 Definitions 
READ page K91 of Section 3-2. 


Notes 


(i) line 8 (equation (3-5)), page K91. This equation can also be written L(y) = h, 
a form which emphasizes the fact that both y and /rare functions and that L is a 
linear transformation from one function space to another. 


Domain C^(1) Codomain C(t) 
(i) fine 10. “ defined on 7" means that the domain of L is C"(/). 


“h is identically zero on Z” means that A(x) = 0 for all x in 7, or in other words 
h is the zero function with domain J. 

(ii) ines 11-12. Normal means that the image of / under the function a, does not 
include 0; i.e. a,(x) # 0, for all x€7. The word “vanishes " is often used to mean 
“is equal to 0", 

(iv) line 13. The reason for requiring y e C'(/) is to be sure that it is in the domain 
of L. “ Satisfies the equation identically on 7” means that Ly(x) = h(x) holds for 
all xe /, or, in other words, that the functions y and 4 (with domain 7) satisfy the 
equation L(y) = A. 

(v) line — 1. “analysis” means the branch of mathematics that derives from the 
concept of a limit; it includes, for example, the differential and integral calculus 
(treated rigorously) and the theory of differential equations. 


19 


Exercises 


1. 


For each of the following differential equations, state whether or not 
it is linear. 

(i) 5Sx*p" + 2y = 73? 

Gi) (ey) = Sxy 
(ii) Q^y + 52 =7 
(iv) y + Sxy? = 0 

(v) Lz 4x (given that y(x) # 0 for all x). 

y 


(vi) The population growth equation y' = ay where a is some real 


number. 
2. For each of the following linear differential equations, state whether 
or not it is (a) homogeneous, (b) normal. 
G) (2e*D2 + xD + e”y = 0 (xe[-1, 1] 
(ii) (2 cos xD? + D)y = 5x (xe[~1, 1) 
(ii) (2cosxD? + Dy =0 (xe[-2,2]) 
(iv) œ- lx" +y =x- |x| (xe[-1,1) 
() G-IxD»" +y =x- ]x] | Gett2) 
3. Which of the following are solutions of the differential equation 
(D? — 4)y =4e* (xe P”? 
() y= 
(i) y = xe 
(Hi) y = (1 + x)e2* 
Solutions 
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l. (ü Linear: it can be written (5x? D? + 2)y = 7x?/2, 
In terms of the Equations (3-5) and (3-6) on page K91, 


L=5x?D? +2 
A(x) = 7x? 
n=2 
a(x) = 5x? 
a, (x) = 0 
a(x) = 2. 


(ii) Linear: it can be written (e*D? + 2e* D + e* — 5x)y = 0. 

(iii) Not linear: the term (y")? spoils it. 

(iv) Not linear: the y? spoils it. 

(v) Linear: as y(x) is never zero, it can be written as 
y’ = 4xy, i.e. (D — 4x)y = 0. 

(vi) Linear: it can be written (D — a) = 0. 


2. (i) (a) Homogeneous, because the A(x) of Equation (3-5) 
is zero. 
(b) Normal: 2e* z 0, for all xe [-1, 1]. 
Gi) (a) Not homogeneous: A(x) = 5x. 
(b) Normal 
Gii) (a) Homogeneous 


(b) Not normal: cos x = 0 when x = t 5. 


Gv) (a) Not homogeneous: x — |x| #0, x <0. 
(b) Not normal: x — |x| =0, x > 0. 
(v) (a) Homogeneous 
(b) Normal: since x — |x| is zero throughout. the 
domain, rather than merely at certain points in it, 


the equation reduces in this domain to a first-order 
normal linear differential equation, namely 


yQ)20 (xe [,2D 


3. (ü) NO: (D? — 4)e2* = 4e?* — 4g?» 
=0 
# 4e, 

(ü) YES: (D? ~4)xe** = D(e!* + 2xe?*) — 4xe?* 
= 2e?* + 2e?* + 4xe?* — 4xe?* 
= 4e, 

(i) YES: (D? — 4)(1 + 3e? = (D? — e? 

+ (D? — 4)xe?* 
by linearity, 

= 0 + 4e?* 

= 4e, 


4.2.2 Homogeneous Equations 


As a start in exploiting the linear algebra covered in earlier units to solve 
differential equations, we consider the special case in which A = 0. 


READ the first paragraph on page K92 (down to “. . . choice of termino- 
logy.”’), followed by Example 1 on pages K92-3. 


Notes 


(i) line 5, page K92. “ of the x-axis"- this phrase merely indicates that the variable 
x will be used later on for elements in the domain. 

(i) line 6, page K92. “Null space" is the term used in K for what we have been 
calling the kernel of a linear transformation: the inverse image of the zero vector 
in the codomain. We discussed the kernel in the Foundation Course and also in 
Unit 2, Linear Transformations (see page N31); the discussion in K (with some new 
examples) begins on page K55, but you need not study it unless you wish to 
revise and have time to do so. 

(iii) line 5, page K93. The solution set defined by Equation (3-11) can also be 
written sin x, cos x) (see page N12). 


Exercise 
(i) How many vectors are needed to span the solution space of 
(D? — 4)y = 0? (Domain of y is R.) 
(ii) Show that the functions y = e?* and y = e~?* (with domain R) are 
solutions of the differential equation in (i). Are they linearly inde- 


pendent? 
(iii) Write down the general solution of the differential equation in (i). 


Solution 

(i) Since the equation is normal, homogeneous and has order 2 
(the operator D? — 4 contains D? but no higher power of 
D), the solution space has dimension 2. 

(i) They are solutions, because (D? — 4)e?* = 0 and (D? — 4) 
e^?* = 0. They are also linearly independent. To prove this 
we show that 

ce c0 7*0 — (xeR) 
implies c; = c3 = 0 (compare the first 3 lines of page K93). 
If, for example, we set x = 0 and x = | in the equation, then 
ce, +c; = 0 
ce? 46e? = 0 


which, solving for c, and cz, give e, = c; = 0. 
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(ii) Since e?* and e~?* are linearly independent and the solution 
space has dimension 2, these functions form a basis. Thus 
the solution set is ¢e?*,e~?*> and the general solution is 
y =c,e?* + c; e7?*, where c, and c, are arbitrary. 


4.2.3 Nonhomogeneous* Equations 
We saw in Unit 3, Hermite Normal Form (page N63) and also in various 
places in the Foundation Course that the solution set of any linear problem 
Ly=h 
is given by 
solution set = {po} + (kernel of L}; 


where yg is any solution of the given linear problem, and the notation 
means that the solution set consists of all vectors (in the domain of L) 
having the form 


Yo + an element of the kernel. 


The next reading passage explains how we can exploit this particularly 
simple structure of the solution set. 


READ the remaining paragraph on page K92 and Example 2 on page K93. 


Example 

Find the particular solution of 
y+ ypax 

such that y(0) = 2, y'(0) = 0. 

From the reading passage, any solution has the form 
J(x) = x + c, sin x + c; cos x 

Cı, C; arbitrary. 

So that 
(x)= 1 + c, cos x — c; sin x 

The two given conditions determine c, and c3. Thus 
z0) — c; =2 
y(021-4e207.c62-1. 

Thus the particular solution we require is 


y(x) = x — sin x 4- 2cos x. 


Exercises 


l. Write down the general solution of the differential equation 
(D? — 4)y = 4e?*. You will find useful information in the answers to 
Exercise 3 of sub-section 4.2.1 and the Exercise of sub-section 4.2.2. 

2. Find the particular solution of (D? — 4)y = 4e?* with the properties 
x0) = 1, y'(0) = 3. 


Solutions 


1. By Exercise 3 of sub-section 4.2.1, the equation (D? — 4} = 
4e^* has y= xe2“ and y = (] + X)e?* as particular solutions, 
and by the Exercise of sub-section 4.2.2, the associated homo- 


geneous equation (D? — 4)y = 0 has y = c,e?* + C5 € ^?" as its 


* We follow K in dropping the hyphen. 
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general solution. The general solution of (D? — 4)y = 4e?" is 
therefore 


y = xe2* + ce? + c, e7?* 
or 
y — (1 + xe + ce? + c, e7?7. 


The two general solutions are equivalent: the second can be 
written 


y = xe? + (1 + c)e? + c, e7?* 


which has the same form as the first since it does not matter 
whether we call the arbitrary constant c, or 1 + c. 


2. We want to choose the constants c,, c; in the general solution 
Y) = xe? bc, ce?" — (xeR) 
to satisfy 
X0) = 1, y'(0) = 3. 


Now 

X0) = 0 + c + c; 
so that 

Cte, =1 
Also 


y' (x) = (1 + 2x)e?* + 2c,e?* — 2c, e7?* 
therefore 


y'(0) = 1 + 2c, — 2e; 


ce, — c; = 1, 


from which c, = 1 and c; = 0; so that the required solution 
is y(x) = (1 + x)e2* 


If you feel you need more practice on this type of problem, try any of 
Exercises 3, 4 and 7 on page K95. (In Exercise (7a), sinh x means 
3(e* — e^") and cosh x means }(e* + e7*).) 


4.2.4 Nonlinear Equations (Optional) 


The discussion from “ Before leaving this section. - - " on page K93 to the 
end of the section on page K94 shows what goes wrong if we try to apply 
the theory to a nonlinear equation. Similar equations were discussed in 
the Foundation Course (Unit M100 24, Differential Equations I). They are 
not essential to this course, so you need not study the passage. 


4.2.55 Summary of Section 4.2 


In this section we defined the terms 


linear differential equation (page K91) 


homogeneous* (page K91) 
nonhomogeneous* (page K91) 
normal (page K91) 
solution* (page K91) 
null space (kernel) (page K55) 
general solution* (page K92) 
particular solution* (page K92) 
Theorems 


1. (Page K92) 
The dimension of the null space of L, a linear differential operator, 
equals the order of L. 


2. (Page K92) 
The solution set of a linear differential equation Ly = h is 


{y: Ly = h) = {yp} + (y: Ly = 0} 
where Ly, — h. 


* These terms we introduced in Unit 3, Hermite Normal Form. 
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4.3 FIRST-ORDER LINEAR EQUATIONS 
4.3.0 Introduction 


The rest of this unit deals with the first-order linear differential equations, 
one of the simplest types of linear differential equation. The main ideas 
have already been discussed in the Foundation Course (Unit M100 24, 
Differential Equations I); in this section we will revise these ideas and also 


learn about some of the manipulations that are necessary to solve such 
equations. 


4.3.1 Homogeneous Equations 


As in the general theory studied in Section 2, we begin with homogeneous 
equations of the form 


Ly=0 
in which we restrict L to be first order; that is, it has the form 
L = a,(x)D + ao(x) 


If the equation is normal, then since it is first order, its solution space is 
one-dimensional and so it has a basis consisting of a single function; we 
can solve the equation completely if we can find such a basis, i.e. a single 
solution. The first reading passage shows how to find such a solution, 
using the method of separation of variables which you met in the Founda- 
tion Course (Unit M100 24, Differential Equations I), and how to obtain 
the general solution from it. 


READ Section 3-3 from page K95 to “. . . where c is an arbitrary constant” 
in the middle of page K96. 


Notes 


(i) line 6, page K96. Why do we need a,(x) # 0? Can we be sure that this con- 
dition is satisfied? If a,(x) = 0 somewhere in J, then we would be dividing by 
zero for that value of x in the rearrangement of Equation (3-17); we avoid this 
possibility by only considering normal equations, which are just those equations 
with leading coefficient a,(x) non-zero everywhere in the interval J. 


(ii) line 7, page K96. When writing the symbol ; we make the assumption that 


(x) #0 for all x in J. Consequently the separation of variables will not yield any 
solution that violates this assumption, in particular the solution y =0 (the zero 
function). We shall, however, recover any solution or solutions left out at this 
stage, when we construct the general solution at the end of the calculation: at 
this stage all we need is any one solution. 

(ii) line 9, page K96. In the notation of the Foundation Course, the integral 
on the right-hand side of this equation would be written 


a j ao(x) 
— or |x— — 


a aix) 


or, in “abused ” notation 


ac(x) 
a(x)" 
The notation used in K is the Liebniz notation. An introduction to Leibniz 
notation was given in the Foundation Course (Unit M100 13, Integration II), and 


a summary of it is given in the supplementary handbook TI*. 
(iv) line 9, page K96. On the left-hand side of the equation, the integral 


Í 12 dx has been evaluated using integration by substitution. The principle 
y dx 


behind this technique was explained in the Foundation Course unit referred to 


* TI stands for Techniques of Integration. 
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LM 4.3.0/4.3.1 


above. Here all we need is the formula using the technique, which is reproduced 
in paragraph 6 of the section Rules of Integration of TI. The formula is 


Í SOW (x) dx -f f(u) du 


. 1 
where 4 stands for y(x). Taking / to be the reciprocal function x -——> —, we 


x 
l o. _fl xel 
[ea = [i9 ER 


The integral at the right is a standard form: since 


£ 
dx 


get 


(in x) = 1 
x 
we have 
Í L du=Inu (u € suitable domain) 


But what is the suitable domain ? Since In z is defined only for positive u, we can 
take the domain to be at most R+. Thus we obtain, putting back y(x) for u 


j x3 yG)dx-2iny&) OG) ER). 


If we assume, therefore, that y(x) is positive for all x € J, then the left-hand side 
of the equation on line 9 can be taken as In y. We suggest that you amend the left- 
hand side of the equation to In y and that of the next to y. The modulus signs are 
not necessary (they are put in to allow for solutions where y(x) < 0 for all x € J, 
but it is not necessary to consider them) since a positive solution does exist and 
we only need one solution. It is important to realize that we need just one solu- 
tion at this stage: see next note. 

(v) line 11, page K96. The theorem referred to is: the solution space of a normal 
homogeneous linear differential equation of order n has dimension n, 

Here n is 1 and so the single solution we have found is a basis for the entire solu- 
tion space; this space is in fact 


CC f289)). 


The general description is often difficult to follow, so we give an example of a 
particular case. 


Example 
Find the general solution of 
y + ay = 0, 
where a e R and the domain of y is R. (ye C'(R).) 


This is of the form of Equation (3-17) on page K96 with a)(x) = 1, ao(x) = 
a, and so the solution is given by the formula at the end of this reading 
passage. The integral in that formula yields 


G) dx - fax = ax. 


(We choose zero as the constant of integration here, because we only 
need one solution of the equation as a basis for the solution space. The 
arbitrary constant will go in later when we express the entire solution 
space in terms of this basis.) 


Substituting the integral we have just evaluated into the formula at the 
end of the reading passage, we obtain 


YQ) = ce" > 


as the general solution; i.e. the solution set is <e~**). For a = 1, the solu- 
tion set is illustrated in the diagram. 
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At the end of any calculation, suchas the one we have just done to solve the 
equation )" + ay — 0, it is advisable to check the result. Here we have 


y'(x) = —ace^** 
ay(x) = ace" "* 


so that y'(x) + ay(x) = 0 and the result is verified. We suggest that you 
work through this example following the steps from Equation (3-17) to the 
end explicitly. 


Exercises 


|. Find for dx, by using the substitution u = x?. 


2. Exercise 1, page K99. Take the domain of y to be R*. 
3. Exercise 3, page K99. Take the domain of y to be (0, 4x], and use the 
substitution u = sin x. 


Solutions 


1. The formula for the substitution is 


IL ((x))y (x) dx = frw du 
where u = y(x) which, in this case, is x?. Therefore 

y(x) = 2x. 
If we now take f = exp, then we have 

feroa -2x dx = [ew udu 

= exp u. 

The only difference between this and the integral we want to 
evaluate is the 2 in the integral on the left-hand side; so we 
conclude 


|: exp (x?) dx = 1 exp u = 1 exp x°. 


It is worth checking the answer by differentiation. 
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xy(x) 42x) =0 (xe R*) 


Here a,(x) = x, ao(x) = 2, using the notation of the reading 
passage. The integral we need is 


[Ea = |Ë. 
a,(x) x 


since x e R*. 


The general solution is therefore 
als c + 
X(x) = ce == CGem') 


(Remember exp and In are inverse functions; In has R* for 
domain.) 
(sin x)y'(x) + (cos x)y(x)=0 (xe (0, x) 


The integral we need is 


ao(x) _ [cosx 
| a,(x) dx= un. = (ee Oda). 


The suggested substitution is u = sin x. To use this we apply 
the substitution formula. The formula becomes 


In x) cos x dx = IL du 


and with f the reciprocal function x ———> E we find that 
x 


c 1 
IE: dx- iF du = In u = In (sin x), 


since sin x is positive for all x in the domain (0, 42]. Thus 
the general solution is 


VOLI ein) = x (x e (0, 1x]) 


4.3.2 Nonhomogeneous Equations 


Now that we know how to solve the Beneral homogeneous equation of 
first order (that is, how to reduce it to the evaluation of an integral), we 
can also solve any nonhomogeneous equation for which we can find a 
particular solution. The next reading passage describes a general method 
for obtaining the necessary particular solution. 


READ from “To obtain a particular solution . . ." on page K96 as far as, 
but not including, Example 3 on page K98. 


Notes 


(G) line —8, page K96. The basic idea of the method is the same as for the 
homogeneous equation: we multiply both sides of the equation by a factor which 
makes it possible to integrate both sides even though we do not know the function 
X. Such a factor is called an integrating factor. The factor 1/[a,(x)y(x)] which 
worked for the homogeneous equation no longer works now, because we cannot 
integrate A(x)/[a,(x)y(x)] without knowing y. Fortunately, even for the non- 
homogeneous equation there is a fairly simple integrating factor, 


a(x) 
exp f a) dx, 


for the equation we are studying. In K this integrating factor appears from no- 
where, but there is a plausible argument for obtaining it; since this argument is 
given in Unit M100 24, Differential Equations I, we shall not treat it in this 
course. 

(i) line 5, page K97. This is an important point: it is much better to understand 
how a method works and be able to work out or look up the formulas when they 
are needed, than to try to memorize complicated expressions like Equation (3-20). 
ue the main thing to remember is that an integrating factor for an equation of 
the form 


YQ) + Py) = Q(x) 


a(x) 


is exp ILS dx. In using the procedure you must always remember to divide 


through by a;(x) first. This is always possible if the equation is normal. It is not 
even necessary to solve the homogeneous equation first, as the examples that 
follow show. 

(ii) line —11, page K97. Notice the difference from the general method 
described at the beginning of the reading passage. Instead of completely solving 
the associated homogeneous equation first and then finding a particular solution 
of the nonhomogeneous equation separately, we do both parts at once by in- 
cluding a constant of integration c when we do the main integration (although we 
choose zero as the constant of integration in the integration that gives the 
integrating factor—e*? here—because one integrating factor will do to solve 
the equation). 


(iv) line —9, page K97. For | xe", see Exercise 1 of sub-section 4.3.1. 


(v) line —5, page K97. The domain of y does not have to be (0, co) or (— co, 0); 
it could be any interval that is a subset of one of these. 

(vi) line — 1, page K97. The use of |x| where you might expect x is a trick for 
solving two separate problems with the same calculation: one problem referring 
to the interval (0, oo) (i.e. when x > 0) and one to the interval (— co, 0). (See 
Section III.2.1 of TI.) 

(vii) line JJ, page K98. “ defined on the entire real line" means “with domain 
R”. 

(viii) line 13, page K98. “ potentially misleading ". Note the warning. It is often 
convenient to forget about the interval on which the functions considered are 
defined and whether or not the functions are continuous. This is all very well 
until something goes wrong, but even if nothing does apparently go wrong, we 
might have an error and we cannot be absolutely sure of our answer unless we 
have applied the theory correctly, that is if we have used the theory only in those 
situations where we know it can be applied. 


You may have realized that all the examples in this section are carefully 
contrived to give simple results. In real life, things are not so easy. It 
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may not be possible, for example to do the integrations; thus, in Example 
1 on page K97, if the right-hand side of the differential equation had been | 


instead of x, we would have had the integral fe dx (in place of fe dx), 


and this integral cannot be expressed in terms of elementary functions.* 
Another thing that can give trouble is that, even when the integral can 
be done, the formula it leads to may be too complicated to be useful. In 
either case, it may be necessary to find some other method (perhaps a 
numerical method) for solving the equation. We shall be looking at some 
of these numerical methods later in the course. An obvious modification, 
which will help in some cases, is still to use the method described as far 
as the integral form of the solution in Equation (3-20), and then to use a 
numerical method, such as Simpson's rule, to evaluate the integrals over a 
specified interval. 
Exercises 
l|. Find the general solution of 

2y()-yx)-€ ER). 
2. Find the solution of 

2y(x)-y(xX)2x  (xeR) 


that satisfies y(0) = 0. 


Solutions 
1. 2y'(x) + y(x) = e* (xe R) 
Divide by a,(x): 


a 1 l x 
ya) +; yQ) = ex. 
An integrating factor is 


1 
exp E dx — exp 6 x) . 


The equation becomes 


yed 4. 5 reset = ; gn. 


d 1 
og 1/2 — ` 9(3/2)x 
p Oye ) >° 3 
Integration gives 


1 
yG9get ms = 2 Jens dx +c 


yeI = — gx c 


wie 


Hence the solution is 
" 1 x -(1/2)x 
yO) 73e + ce (xe R). 


* The elementary functions are x", In x, exp x, sin x, cos x, etc. 
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2/690 + y(x) = x 
Since the linear operator, 2D + 1, is the same as in Exercise 1, 


the integrating factor and the kernel are also the same as there. 
The method of Solution 1 gives 


d 
de 09e = (ven +3 yo) anm 


so that 


Y(x)et/2)* = ipee dx + c. 


A 1 
To integrate 3 xe/2)*, we use integration by parts. 


[^ Dee fra fex nr 


or in Leibniz notation 


| SO) dx =FR — | O dx. 


Jo) =x 


£(9) = etx 


Formula yields 
d 1 
x = (e25) dx = xes | 2s qx. 
dx 
Therefore 
P xe/2* dy = xex _ Qelt/2)x — (x— 2)et/2«. 


Thus 
yGgett/? = (x — 2)et/2s 4 e 
and the general solution is 
Wx) = x — 2 + ce (2, 
We want (0) = 0, 
ie. 
02 -2-c. 
So the required solution is 
))=x—2+2e 02 3 (xe R). 


If you require further practice, try Exercises 5 and 7 on page K99. 
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4.3.3 Bernoulli’s Equation 


The last reading passage in Section 3-3 of K deals with a special type of 
non-linear first-order differential equation which can be solved by means 
of a substitution converting it to a linear equation. It is less important than 
the rest of Section 3-3, but this type of non-linear equation crops up occa- 
sionally. The method for solving it is the kind of thing you should know 
where to look up, but should not try to memorize. 


READ Example 3 on pages K98-99. 


Notes 


G) line —6, page K98. If n > 0, this step is only valid if y(x) # 0 for all x in the 
domain of y. This is the reason for the statement on line 2 of page K99, that the 
solution y = 0 has been “suppressed ". The procedure may also “ suppress ” other 
solutions of the equation, as the example discussed on pages K93-94 shows 
(that example is a particular case of Bernoulli's equation). 

(i) line 17, p we K99. The choice of c is restricted by the condition that u = y^! 
must not take the value 0 for any x in the domain of y. 


Exercise 
Find the general solution of 
y -Py--y 
where P is a positive number and y: Rt ———> R*. 


(This solution is used in the next passage from K.) 


Solution 


y — Py = —y! is a Bernoulli equation with n = 2. Following the 
prescription on page K98 we multiply by y^? and then write the 
equation in terms of u = y^. Both steps are fully justified since 
the codomain of y does not include 0. The first of these steps gives 
yy — Py = -1 
which is equivalent to — w’ — Pu = — 1, that is, to 
u + Pu = 1. 


An integrating factor for this equation is eP*; this gives 

d 

= (u(x)eP*) = (u'(x) + Pu(x))e?* = e”, 
Integration gives 

1 

u(x)eP* = — ePx 

(x)e?* = pe tc, 
so that 


1 
— a 
u(x) P + ce 


and (by the definition of u) 


1 P 
(I/P) + ce 1 + cPe~?* 


xx) = (xeR*). 


To ensure that the codomain of y is R* as required, we must 
require c > —1/P. 
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4.3.4 Summary of Section 4.3 


In this section we defined the terms 


integrating factor (page K97) 
Bernoulli’s equation (page K98) 


Techniques 


1. We have seen how to solve the following types of differential equation. 


(i) Linear homogeneous and normal 
(a,()D + ao(x))y = 0 
(i) Linear nonhomogeneous and normal 
(a,(x)D + ap(x)y) = h 
(ii) Bernoulli’s equation 


(a,() D + ag(x))y = hy" 
2. We have also seen how to perform 


fres, Je. a#0, 


These results are recorded in TI. 
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4.4 SOME APPLICATIONS 
4.4.0 Introduction 


In this section we consider some examples illustrating how differential 
equations can be used to model real situations, and their solutions used 
to make quantitative statements or predictions about.these situations. In 
these applications, solving the differential equation is only one part of the 
mathematician’s work; he must also set up the differential equation to 
represent the real system he is studying, and interpret its solutions in terms 
of the behaviour of the real system. 


[rar sstem | modelling differential 
real system ——ə—.— I 
I equation 


| 


1 
| solving 


interpretation 
behaviour 


We shall discuss other more complex modelling situations in Unit 9, 
Differential Equations II, Unit 11, Differential Equations III, and. Unit 13, 
Systems of Differential Equations. 


4.4.4 Growth and Decay 


The application of first-order differential equations to population growth 
problems was discussed in the Foundation Course (Unit M100 24, Dif- 
ferential Equations I). The next reading passage covers roughly the same 
ground. 


READ Section 4-9 on page K166, to Figure 4-4 on page K167. 


Notes 


(O) I line 15, page K166. We are making the approximation that the number of 
individuals, a positive integer, can be adequately approximated by a con- 
tinuously differentiable function y. 


Gi) line —1, page K166. The notation here is not quite the same as in the solu- 


tion in the previous section of this correspondence text: c here den 
: otes the num- 
ber denoted by cP there. s 


(ii) line 4 and Figure 4-4, page K167. From these mathematical results we can 
make predictions about the behaviour of any population that satisfies the dif- 
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ferential equation. The most obvious is that when z is very large, e^ ”* is very small 
and so y is very close to P; in other words 


lim y(t) =P. 


t large 


Thus the population approaches the size that can just be supported by the 
available supply of the necessities of life. 


Exercises 


]. Exercise 9, page K169 (express your answer in terms of the relevant 
constant of proportionality). 
2. Exercise 4, page K169. 


Solutions 


l. Step I: setting up the equation 
Let the raindrop have radius r(f) at time t. Then its volume is 


4 
+ [r(0)]° and its surface area is 4x[r(1)]?. 
Evaporation causes a rate of decrease of volume proportional 
to its surface area, i.e. 
d (4x a à 
SSF ror) = -k -sron 


where k >0 is a constant of proportionality. The equation 
can be written 


4n(r(OPr'(t) = —k4x[r(0Ë 
or 

r'(t) = —k, 
provided r(t) # 0. 
Step 2: solving the equation 


It can be integrated at once (i.e. the integrating factor is 1) 
which gives . 


r(t) = —kt + c. 
Step 3: interpreting the solution 
We are told that 

r(0) = ro, 


OR c-r 
and are asked to find t, such that 
r() =0 


—kt, + c = 0. 


The first equation tells us that c = ro and the second then tells 
us that t, = ro/K. 


2. Step 1: setting up the equations 
The equation of growth is given: 
y= (ki — Kk». 
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We must also consider the equation that would hold if there 
were no births, i.e. if k, = 0; this equation is 


y = —k; y. 
Step 2: solving the equations 
The general solution of the growth equation is 

y(t) = cet 0 
That of the no-births equation is 

X) = de^ ^it (2) 
Step 3: interpreting the solution 


We are told that the colony doubles in size every 24 hours; 
therefore, with t measured in hours, we must have (by Equa- 
tion (1)) 


ce7*2^ 2 units of population 


ce*17*39 — ] unit of population 


therefore 


gi 24 — 2 


k,—k;-4;In2 


We are also told that if there were no births, the colony 
would halve in 8 hours, so that as above (by Equation (2)) 


e kes — 1 


k, = 41n2 (3) 
To find k, we substitute for k; in Equation (3), obtaining 


k, = š In 2 = 0.12 (correct to two decimal places) 
k, = $ In 2 = 0.09 (correct to two decimal places) 


that is, the growth rate is about 12% per hour and the death 
rate is about 9% per hour. 


4.4.2 Electric Circuits 


The formulation of differential equations for electrical circuits containing 
resistances and condensers is treated in the television programmes of this 
and the previous unit, A summary of the rules is given at the top of page 
K170, which includes also the type of circuit element called an inductance. 
The next reading passage discusses a simple circuit containing a resistance 
and an inductance. The physics is different from the circuit considered in 
the television programme, but the mathematics is almost identical, as the 
following diagrams and equations show. 


T R 
v(t) — L 
vc sotage R battery 
TV K 
R eam VO) LS Ri= V() 


Equivalent roles are played by the following pairs of terms: 


charge on condenser g : current i 
resistance R: inductance L 


; 1 š 
inverse capacitance T resistance R 


READ Section 4-10 on page K170 to‘... term "impedance" here” in the 
middle of page K172. 


Notes 


(i) law (b), page K170. The “voltage drop” across an electrical component is the 
difference in electrical potential, as measured with a voltmeter. To law (b) we 
might add the rider that if a voltage source is applied to an open loop of the circuit 
then the sum of the voltages across the various components is equal to the applied 
voltage: 


V= V. + V. 


S 
Applied - 
Voltage e 
= 
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(i) line —3, page K170. The “easy computation" is a particular case of 
Exercise 7 on page K99. A differential equation with a condition of the type 


i(0) = something 


is known as an initial value problem. These will be discussed in detail in Unit 9, 
Differential Equations II. 

(iii) line 8, page K171. The “alternating current source" has a time-varying 
voltage of the form Esin ot (E constant). The differential equation is to be 
solved on the interval J = [0, co). 

(v) line 12, page K171. The calculation leading to Equation (4-92) is given as an 


exercise below. . 
(v) lines —9, —8, page K171. The definitions of Z and « are shown in the dia- 


gram. Notice the use of the formula 
sin (0 + é) = sin P cos é + cos 0 sin ¢, 


with 0 = wt and ó = —«. 


wl 


Exercises 


1. Solve Equation (4-91), page K171 (take 1 e RZ). 
Exercise 3, page K175 (the voltage across a condenser is g/C where q is 
the charge, related to the current by i = dq/dt). 


Ignore the coulombs, microfarads, seconds, and ohms, but take 
C = 300 x 1075, not 300. 


3. Set up the differential equation for the following circuit: 


Alternating current source 


J 


Solutions 


1. Applying the method of Section 3-3 of K (pages K96-7) we 
tewrite Equation (4-91) as 


di zi R, E diii 

S++ i= ssi 2 

di p cpm 
An integrating factor is 


exp Jaw dt = exp (R/L)t; 
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using it, we obtain 
G + R i) ex Ree |, Bi 
dL po oa (eet) 


= (sin wt) Rt 
=F exp T: 


Integrating both sides and using Section IIL5.2 of TI we 
obtain 


LN 3 

R. L PT 7 sin at — o cos or) 

rp =? SE 
pte 


so that the general solution is 


E/R, 
m Z sin ot — o cos wt) Ri 
R i + cexp/ i). 
pte 


To pick out the solution required, we determine c from the 
condition i(0) = 0. The general solution gives 


LE 
L 
0=i0) = te 
(e +=) 
so that we obtain 
E [^] 
L 


substituting this into the general solution, we obtain 


L 


i= R + 2 
pte pte 


z (finos w cos wt) 7 wexp( *) 
] p^ 


t which is equivalent to the corrected version of Equation 
(4-92). 


(a) Step 1: setting up the equation 
Kirchoff's law (b) gives 
iR + s 0 
i ci^ 


and since i — 2 the differential equation of the circuit 


(with the switch closed) is 
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Step 2: solving the equation 


Separation of variables (see page K96) gives the general 
solution 


Step-3: interpreting the solution 


We are given g = 0.03 at t = 0 and so c = 0.03 and the 
required solution is 


—t -t 
q = 0.03 exp Rem 0.03 exp reg 
(b) 


q(t) 
0.03 
0.02 


0.01 


Alternating current source 


oa 


Kirchoff’s law applied to point A on the diagram gives 


A(t) + i(t) = io a) 
Kirchoff’s second law applied to closed loop ABC gives 

Ri, (t) — Cq(t) = 0 (2) 
where 

; = dq; 

i) = T (3) 


We must reduce Equations (1) and (2) to a form which we 
can solve. 


Differentiate Equation (2): 


so, using Equation (3) 


i y 
j^ Cialt) = 0. 


R 


Substitute into Equation (1) for i(t): 


di 
R a + Ci,(t) = Cis. 


Please note that for examination and assessment purposes you are not 
expected to be able to recall Kirchoff's laws and the form of voltage drops 
across network components. But, given the laws and the form of the voltage 
drops, you should be able to set up the appropriate differential equations. 


4.4.3 Summary of Section 4.4 


To analyse a physical situation in terms of a differential (page C35) 
equation we follow the procedure: 


step 1: set up the equation 
step 2: solve the equation 
step 3: interpret the solution 


Kirchoff's Laws enable us to set up the differential equa- (page K170) 
tions of an electrical network: 


Law (a): The algebraic sum of the currents flowing into 
any point in an electrical network is zero. 

Law (b): The algebraic sum of the voltage drops across 
the various electrical components in an oriented 
closed loop of a network is zero. 


The voltage drops across the electrical components con- (page K170) 
sidered in this section are given by 

g/C for a condenser of capacitance C 

iR for a resistance R 
di 


Un 


for an inductance L 


where g is the charge on the condenser and i is the current 
through the resistance or inductance. 


The rate df accumulation of charge on a condenser is (page K172) 
given by 
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4.55 SUMMARY OF THE UNIT 


We have looked at a particular infinite-dimensional vector space—the 
space of all real functions—and have chosen suitable subspaces of ade- 
quately differentiable functions. On these subspaces we have defined 
linear transformations of the form 


a,(x)D" + a, ,(X) DD"! +  ag(x). 


Such linear transformations are called linear differential operators, and a 
linear problem involving such an operator is usually referred to as a linear 
differential equation. In the special case when the functions a,(x), ..., a(x) 
are constant functions, we have a constant-coefficient linear differential 
operator. (Their use will become apparent in Unit 9, Differential Equations 
I) 


The solution space of a linear homogeneous differential equation is finite- 
dimensional and its dimension equals the order of the equation. The 
terminology for linear differential equations is mostly coincident with that 
of the finite-dimensional linear problem which is given in Unit 3, Hermite 
Normal Form. 


The normal first-order equation, which was solved in Unit M100 24 is 
revised here. Applications which were considered include growth and decay 
of populations and carefully chosen electric circuits. 


Definitions 


The terms defined in this unit and page references to their definitions are 
given below. 


continuously differentiable ^ (page K13) 


interval (page K3) 
linear differential operator (page K86) 
order (page K86) 
linear differential equation (page K91) 
homogeneous (page K91) 
nonhomogeneous (page K91) 
normal (page K91) 
solution (page K91) 
null space (kernel) (page K55) 
.general solution (page K92) 
particular solution (page K92) 
integrating factor (page K97) 
Bernoulli's equation (page K98) 

Theorems 

l. (Page K92) 


The dimension of the null space of L, a linear differential operator, 
equals the order of L. 


2. (Page K92) 
The solution set of a linear differential equation Ly = h is 


{y: Ly = h} = {y,} + (y: Ly = 0} 
where Ly, = h. 


Techniques 


1. Use (aD + b)(cD + d) = acD? + (ad + bc)D + bd when a, b, c and d 
are constant functions. 
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2. When a, 5, c and d are not constant functions, apply the operator in E 
1 by applying each factor in turn. 
3. Wehave seen how to solve the following types of differential equation. 


() Linear homogeneous and normal 


(a,(2)D + ao(x))y = 0 . 
Gi) Linear nonhomogeneous and normal 
(2,0)D + aq)» =+ ` 


G) Bernoulli's equation 


(a,(x)D + ao(x))y = hy" 


4. We have also seen how to perform 


fres, Je #0, 
j ü E 
xe™, | —. 
» 


These results are recorded in TL 


5. To analyse a physical situation in terms of a differential (page C35) . 
equation we follow the procedure 


step 1: set up the equation 
step 2: solve the equation 
step 3: interpret the solution 


Notation 

ew (page K86) 
Cla, b] (page K3) 
Dip (page K50) 
e) (page K86) 
(a, b) (page C16) 
(a, co) (page C16) 


(— co, b) (page C16) 
(—00, 0) (page C16) 
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4.6 SELF-ASSESSMENT 
Self-assessment Test 


This Self-assessment Test is designed to help you test quickly your under- 
standing of the unit. It can also be used, together with the summary of the 
unit for revision. The answers to these questions will be found on the next 
non-facing page. We suggest you complete the whole test before looking 
at the answers. 


1. Which of the following are linear differential operators? (y e C?(I) 
in each case) 


(a) y—— y* + 3y! «xy 
(b) y— yt”) + 3y" + xy 
(cy y —— py’ + 2p 
(d) y— y + 2y +3 
2. Find (x D)?f, where f(x) = x? (x e R). 


3. Which of (a), (b), (c), (d) and (e) correctly complete the following 
statement ? 


If F is the derived function of Ge C'(I), (n e Z*), then F necessarily 
belongs to... 


(a) Cu) 
(b) C°) 
() @' (D 
(d) CQ) 
() ern 


4. For each of the following equations, say whether it is 


(a) homogeneous on J 
(b) normal on J 

(c) both 

(d) neither 


G) x?y"(x) — xy'G) + 30) = 0 
with I = (—3, 3) 
Gi) x*y"(x) — xy'(x) + 3309 = e7* 
with I = (3, oo) 
(ii) y"(x) + (sin x)y(x) = x 
with I = (—7, z) 


Gv) y"(x) + Gin x)y(x) = 0 


š n 
with I = (0 5) 
() xl'Q) =sin x, 
with I = (—z, 0)° 


5. For each of the following, say whether it is a solution of the equation: 
(D? + D—2)y(x)=0 (xe R) 


G) ya) =e* 
Gi) yx)-e7 
(iii) yo) = e?* 
(iv) sx) =e 


6. Write down the general solution of the equation in Question 5. 
7. For each of the following, say whether it is a solution of the equation: 
(D? + D — 2)y(x) = —2e7* (xe R) 
(i) yx) = e* 


10. 


11. 


12. 


13. 


14. 


Gi) »9)-e7* 
Gü) y) = e? 
(v) YQ) =e" 


Write down the general solution of the equation in Question 7. 
Find the integrating factor for each of the following equations: 


G) 3y'(x) + 4y(x)=e7* (xe R 
(Gi) 2xy')— y(x) = x? (xe R*) 


Write down the general solution for each of the equations in Ques- 
tion 9. 


What is the null space of each of the operators defined by 


G) L:y— 3y' + 4y 
and (ii) L:y—— 2xy' — y, 


the linear operators in Question 9? 
By the substitution u = y?, reduce 
Jy + xy? — 2x = 0 
to a linear equation. State the linear equation. 


For each of the following, state the order of the sum L, + L;. All 
operators are defined on the interval ( —z, z). 

(i) L,=2xD+3,L,=xD—1 

(ü) L,-e"7D'—D-a2 

L,-—e7D't2D42 

(i) L,2D-4L;-sinxD'—D-al 
Write down a differential equation for the circuit below, introducing 
suitable notation. 


Alternating current source 


C) ig sin wt 


L 
0000000 
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Solutions to Self-assessment Test 


1. 


il, 


12. 


13. 


(b) only. 
Compare (a), (b). (c) and (d) with Equation (3-1) on page K86. 


(xDYà = (xD(x Dx?) = (xD)(x3x?} 
2x3. 3x 
-9x3 


Only (c) and (d) correctly complete the statement. (See sub-section 
1.2 of this correspondence text.) 

G (9 G) (b) (i) (b) 

(ivy) © (v) (b). 

(i) YES (ii) NO Gi) NO (iv) YES. 

(Check by substitution.) 

yY =e + ce? 

€,, €; are arbitrary constants. 

(i) NO (ii) YES (iii) NO (iv) NO 
yme7 ree + cx e ?* 

(i) ELS - e^3 (ii) ef 7 12x dx = x 1⁄2 


(i) Using the integrating factor of 9(i), the equation reduces to 


yet = 1 fore 


=e +c 


y= e * + ce Ho, 


(i) Using the integrating factor of 9 (ii), the equation reduces to 


112 1 
- x7 My 
yx maf dx 


2 
3/2 
x 
= |—dx 
j: 
x)? 
CS 
ie. 
3 
a 
pL ext? 
FES 


The null-space or kernel are those functions that map to the zero 
function under L. 
G) From the solution to question. LO(i), the null space is (e^ 4/3», 


(à) From the solution to question LO(ii), the null space is x'/2y, 
Ius’,  w-2yy 
Hence 

lu + xu — 2x = 0. 


() First order 
(à) First order (L, + E; = B + 4) 
(i) Second, order 


14. 


Alternating current source 


O io Sin wt 


L 
^ 0000000^ -— 
4 
A B 
£ ANAAA^ 
l2 R 


At a node 
lj + i> = io Sin ot 
where iy and o are constants. 


Voltage between A and B is 


LR (via R) 
di, š 
£ x (via L) 
Hence, 
di 
iR=L PE 


So, eliminating i; 


di TNT d 
L a = Rig sin wt — ij) 


di š 
L E + Ri, = Rig sin wt. 


47 APPENDIX (OPTIONAL) 

A differentiable function f, whose derived function is discontinuous yet has 
the same domain as f. 

Let f be the function defined by 


- I 
fxh sin Ix d Eie (x e R) 


Then f has the graph 


This function is differentiable at every point in R and 


2x sin 1/x — cos 1/x if x0 
0 if x20 


But Df is not continuous at 0 : Df(0) = 0, but lim Df(x) does not exist. 
x0 


Dfix— | 


First, we show that Df(0) = 0. We go back to first principles: 


£0 -s0 
h-0 


ENDE Uh — 0 
h^ h 


Df(0) = lim 


= lim A sin l/h 
h^0 


Given e > 0, consider ñ e [—6, £] ^ #0. 


Then, 
T patx|sintl<yat i intl ct 
l|- in-|« R is 
; h sing since sin — | < 
ie. 
E 
hsin -e [—e, €]. 
h 
Therefore 
š ek 
lim hsin- 20 
hoo h 
so that 


Df(0) = 0. 


Next, we show that Df — f' is not continuous at 0. 
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For f” to be continuous at x = 0, we would have to be able to choose, for 
any £ > 0, an interval [— ó, ó] such that whenever x is in this interval, 
If’) —F'O)| < s; that is 


I| < e. 
But 
r 1 cw 
f'(x) = —cos Ë + 2x sin É (xe R, x #0) 
x x 
and although the second term gets as near to zero as we please for, x close 
en 1 " 
to zero, this is not true of the —cos () term, which oscillates between 
x 


— l and 1, as x gets closer to 0. So it is intuitively clear that this function 
is not continuous at 0. To be rigorous about it, we choose £ = 4 (any value 
less than 1 would in fact do), and proceed to show that we then cannot 
choose an interval [— 6, 6] with the required property. So let 5 be any 
number greater than zero, and let us examine the interval [— ó, ó]. If we 


: Š š 1 
choose an integer N with Nz > 1/6, then letting x= we see that 
T 


x € [— ó, ó], and 


IA I z f 
f(x) = —cos (3) + 2x sin (3) 
= +1 +0. 


Therefore 
Ife) 21»:-t 


and f” is not continuous at 0. 
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